Model-Independent Predictions for Low Energy Isoscalar Heavy Baryon Observables 
in the Combined Heavy Quark and Large N c Expansion. 



Z. Aziza Baccouche, Chi-Keung Chow, Thomas D. Cohen and Boris A. Gelman 
Department of Physics, University of Maryland, College Park, MD 20742-4111- 

Model-independent predictions for excitation energies, semileptonic form factors and electro- 
magnetic decay rates of isoscalar heavy baryons and their low energy excited states are discussed in 
terms of the combined heavy quark and large N c expansion. At leading order, the observables are 
completely determined in terms of the known excitation energy of the first excited state of A c . At 
next-to-leading order in the combined expansion all heavy baryon observables can be expressed in a 
model-independent way in terms of two experimentally measurable quantities. We list predictions 
at leading and next-to-leading order. 
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A number of nonperturbative expansions have been used to study low-energy hadronic observables. The heavy 
[ quark expansion has been very successful in describing hadrons containing one heavy quark. In the heavy quark limit 
the heavy quark spin-flavor decouples from the dynamics and hadronic spectrum exhibits an approximate heavy quark 
spin- flavor symmetry. Heavy quark effective theory (HQET) is a consistent framework to study the 1 /rriQ corrections 
to the heavy quark limit Another useful limit is the large N c limit where the number of colors N c of the QCD 

gauge group is taken to be large j7|||. In the large N c limit the baryon sector of QCD exhibits an approximate 
light quark spin-flavor symmetry; eg. as N c goes to infinity the nucleon and A become degenerate |p|-|l3| . For heavy 
baryons, i.e. baryons containing one (c or b) heavy quark it is useful to consider a double limit and a combined heavy 
J> j quark and large iV c expansion In Ref. |I5j-|l7| isoscalar heavy baryons have been considered near the combined 
\£> '. limit. 

The combined heavy quark and large iV c expansion is formulated in terms of the counting parameter A ~ 
h-had/mQ, 1/N C with AhadNc/rriQ arbitrary (here A^ad is a typical hadronic scale). As A approaches zero, the heavy 
baryon spectrum exhibits an approximate dynamical symmetry — a contracted 0(8) symmetry — which connects or- 
bitally excited states of the heavy baryons to the ground state |l4],[l5|. An effective theory describing low-energy 
f^i ■ excitations of heavy baryons can be developed using the counting rules consistent with a contracted 0(8) symmetry 
""p^ \ [pTsfl . The purpose of the present paper is to use the formalism of Ref. jnj to derive model- independent relations 
among observables characterizing isoscalar heavy baryons and their decays. At present, experimental data is available 
only for the excitation energies of the doublet of the first excited state of A c . At leading order in the combined 
expansion this energy completely determines the excitation energy of the low-energy orbitally excited state of Ab as 
well as the dominant semileptonic form factors and electromagnetic decay rates. Such predictions are necessarily 
rather crude as our expansion is in powers of A 1 / 2 . At next-to- leading order an additional observable is required 
to make further model independent predictions. As data for more heavy baryon observables (such as the excitation 
energy of Af, baryon, semileptonic and radiative decay rates of A c and Af, baryons and their excited states) become 
available, the next-to-leading order predictions can be tested against experimental results. In this paper we discuss 
which observables are most sensitive to the NLO corrections, and, thus, can be used to extract two model-independent 
constants which arise at leading and next-to-leading order in the effective expansion. Using the formalism of Ref. |l7| ] 
we derive model independent relations among a number of observables characterizing isoscalar heavy baryons and 
their decays. 

The effective degrees of freedom near the combined limit are the collective motion of the brown muck — the light 
quarks and gluons — relative to the heavy quark. As shown in Ref. p6[ , the energy of the collective excitations is 
of order A 1 / 2 while the energy of the intrinsic excitations of the brown muck itself are of order A . The effective 
operators that excite the collective degrees of freedom can be defined from the corresponding QCD operators in the 
combined limit. The construction of these operators from QCD in a model-independent way was discussed in detail 
in Ref. They are the total momentum P (a constant of the motion) and its conjugate position operator X of 

the entire system. These operators commute up to corrections of order A with another pair of the collective conjugate 
operators p and x which are related to the momentum and position of the brown muck relative to heavy quark. 

We have so far focused on the isoscalar heavy baryons, A c and Af,, and their excited states. The effective Hamiltonian 
for these states can be written solely in terms of the collective operators discussed above. As shown in Ref. jllj ], the 
contracted 0(8) symmetry naturally arises in the combined limit if the counting rules of the collective operators are 
given by, 

(x, X,X Q ,X e ) ~ A 1 / 4 , 
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(p, P P Q , P t ) ~ (1) 

in the sense that the typical matrix elements of these operators between the low-lying heavy baryon states scale as A 1 / 4 
and A -1 / 4 . The effective Hamiltonian based on these counting rules up to terms of order A was derived in Ref. Eg]: 



: (m H + m N ) + co + ( 2{m f +mH) + & + |^ 2 ) + + O(X^) , 

II II II ' II ( 2 ) 
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where m^r is the nucleon mass; mn is the mass of the heavy meson containing the same heavy quark as the heavy 
baryon. These masses are measurable quantities while the corresponding brown muck and the heavy quark masses 
are not. In the combined limit, the difference between the heavy quark and the corresponding heavy meson mass is 
of order A ; the difference between the brown muck mass and the nucleon mass is of order A as well. The heavy 
meson and the nucleon masses scale as A -1 in the combined limit (r|. The reduced mass, [Aq, which is of order A -1 , 
is given by, 



m N m H ,_, 

MQ = : ■ (3 

m N + m H 

In eq. (||), 7i\n refers to the piece of the Hamiltonian whose contribution is of order A". The term W^i/a is referred 
to as leading order (LO), and TL\ as next-to-leading order (NLO). The combined expansion is in powers of A 1 / 2 . 

The effective Hamiltonian in eq. ^ contains three phenomenological parameters — Co, k and a which are of order A . 
At order A these constants are independent of the heavy quark flavor. However, they contain l/m<g ~ A corrections. 
As seen from eq. the I/tuq corrections to k and a contribute only at order A 3 / 2 , i.e. at next-to-next- to-leading 
order (NNLO) and higher. Hence they can be neglected when we work only up to NLO. The 1/toq correction to cq 
contributes at order A, i.e. at NLO. However, the dynamics of the collective degrees of freedom does not depend on 
the constant Cq which only determines overall ground state energy. Thus, only two flavor-independent constants up 
to NLO determine the low-energy dynamics of the effective degrees of freedom — the collective motion of the brown 
muck relative to the heavy quark. They can be extracted from data in a model-independent way provided at least 
two relevant observables are measured. 

The phenomenology of the A c and At baryons and their first orbitally excited states based on the effective Hamil- 
tonian in eq. (|J) was discussed in Ref. fL7f . In addition to the Hamiltonian, the effective operators determining 
semileptonic form factors and electromagnetic decay rates have been derived in a model-independent way in the com- 
bined expansion [M. Remarkably up to NLO, these effective operators do not contain additional phenomenological 
parameters. In what follows, we present the model-independent predictions of the effective theory at LO and NLO. 
At LO, these predictions are the same as obtained from the bound state picture of the heavy baryon in which the 
heavy baryon is thought of as the bound state of the heavy meson and an ordinary baryon (eg. the nucleon) [ft8|-p7||. 
However, the formalism of Ref. |rjj ] puts these results in a model-independent framework. 

A useful eigenbasis for the low-lying heavy baryons is spanned by the eigenstates of the leading order terms of the 
effective Hamiltonian, i.e. |Aq; N, I, m; <r), where A, I, m are the quantum numbers of the three-dimensional harmonic 
oscillator and a is the third component of the heavy quark spin (a = ±1/2). The heavy quark spin decouples from the 
dynamics up to NNLO. As a result, the total angular momentum of the heavy baryon is J = I + sq = \l ± 1/2|. The 
heavy baryons are the eigenstates of the total angular momentum, i.e. they are appropriate linear superpositions of 
the states |Aq; N, I, m; a). We denote these states by |Aq; N, J, J z ). The ground states in both the charm and bottom 
sectors, A c and A^, have been observed. In addition, the doublet of the first orbitally excited state — (A c i, A* x ) — has 
been observed |28|-|3C|]. For further convenience we use the following notation for these states: 

|A Q > = |A Q ;0,-, J,) ~A C , A h , 

\K Ql ) = \K Q -l,\, Jz) ~A C (2593), A b (?) , 

|A^) = |A Q ;l,|j z ) ~A C (2625), A b (?) . (4) 

The states of the doublet are degenerate in the combined limit up to NNLO. This leads to an order A ambiguity in 
the determination of the energy of the first excited state. One way to define this energy is to take the spin-averaged 
mass of the doublet: 

12 

TOa q = 3 TOA «! + 3 mA Qi ■ ( 5 ) 
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In the charm sector the mass is approximately 2610 MeV. Similarly, the heavy meson mass, can be taken 
as the spin-averaged mass of the corresponding pseudoscalar and pseudovector mesons: 

1 3 
mp, = -m D + -m D , sa 1980 MeV , 
4 4 

1 3 

m B = -m B + -m B - ^"oilQMeV . (6) 
The excitation energies up to NLO are given by, 

Am c Em i ,-m At = 1 p+^-+0(A 3 / 2 ), (7) 
V A*c 4! Kfi c 

Am b = m- K -m Ab = f± + l-SL + 0(A 3 / 2 ) , (8) 

where the order A corrections are determined by treating the a x 4 /4! term in Hamiltonian perturbatively. As indicated 
above, these excitation energies are completely determined in terms of two constants, n and a. 

In addition to the spectroscopic observables, the effective theory predicts the dominant semileptonic form factors 
and electromagnetic decay rates of the excited heavy baryons. Up to NLO these observables are completely determined 
in terms of k and a. Semileptonic form factors and electromagnetic decays have yet to be measured. 

As discussed in Ref. jlTj , the effective theory in the combined limit can reliably predict the dominant semileptonic 
form factors only for the transfer velocities, \5v\, of order A 3 / 4 . However, this is the range of the velocity transfers 
within which the form factors significantly change from their values at zero recoil. 

In the combined limit the hadronic amplitude of the Af, — > A c £v decay up to NLO is given in terms of a single 
independent form factor |l7f| : 

<A c (^)|c 7 ^(l~75)fo|A 6 («)) =e(z)u c (v>)Y l (l-l 5 )u b (v)(l + 0(\)) , (9) 

where Dirac spinors are normalized by uq(v, s)uq(v, s) = 1. 

The form factor O(z) in eq. (|^) is given as a function of the kinematic variable z (instead of the more usual velocity 
transfer \bv\ or momentum transfer q) defined by: 

m N \Sv\ 



As a function of the velocity transfer \5v\ the form factor O has an essential singularity in the combined limit, while 
O(z) is a smooth function. The new kinematic variable z is easily related to the more familiar variable the velocity 
transfer parameter w = v ■ v' where v and v' are the 4- velocities of the initial and final states: 



m Ny /2(w-l) ^ x<nri 3/2- 



« = 7 rr-\ ^ W2 (1 + 0(A^)) • (11) 



The relation in eq. (|Tl]) is valid for velocity transfers of order A 3 / 4 . 
The form factor O(z) was calculated up to NLO in Ref. p7| ]: 

(^ + W) 3/2 eXP (^J 

(i + w^Mvm ( "t^g + 5 ^ - H ) C 1 + °( A3/2 )) • ( 12 ) 

Useful quantities extractable from experiment are the value of Q(z) at zero recoil and the curvature at zero recoil (the 
slope of O(z) at z — vanishes). Up to NLO, these observables are given by |l7| : 

e ° - e <* = °) = t^#W ( 1 + 5 1fi 3/ 4 2 5( ^-7^V0 i 1 + • ^ 

- = **(z = 0) = - 4^Ct, fl + ^ 45( ^'^ 2 ' 160 ^ f 1 + O(X^)) . (14) 
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At leading order, Oo is independent of k and a; it is an overall normalization which is approximately 0.998. This value 
is very close to the heavy quark effective theory (HQET) value of unity, despite the fact that the present expansion 
is not the pure heavy quark expansion. The curvature at zero recoil and a slope at zero recoil of the form factor O as 
a function of w are related by, 

|% = 1) = -^— p. (15) 

OW JPb + JPc 



Note that p is of order A while d<d/dw(w = 1) is of order A 3 / 2 in the combined limit. 

There are two semileptonic decay channels of A b corresponding to the decays to the doublet of the first excited 
state: Ab — > A c iiv and A b — > A* x £p. As shown in Ref. |17| ], the dominant form factor which determines the hadronic 
amplitudes for both of these channels is given by, 

(A c i(t/)|c7 AI (l - l5 )b\A b (v)) = V3 3(a0tte(tJV(l - J 5 )u b (v) (1 + O(X)) , 

(a: 1 («')|c 7 ^(i-75)&|a 6 («)) = ~(z)u c „(tf)(^ 5 ~gnMv)(i + o(\)) , (16) 

where the Rarita-Schwinger spinors are normalized by U u (v, s)u v (v, s) = — 1. 

For the velocity transfers of order A 3 / 4 the form factor 2(z) is determined in Ref. JlTj up to NLO: 

- y z > - K i /4 (v7 u +v ^ ) 2 cxp ^ 

( i j a ( (105Mb-230yMt,A'e+45^e)K 13 1/2 2 1 (1 j_ (n(\3/2\\ (1 J\ 

y 1 + 4i^75(7^+7^) ^ ib^M^c + 2" K * "4 2 J) + )) ■ v-v 

At zero recoil the form factor S(z) vanishes. The slope of E(z) at zero recoil is given by, 



95 4 Mb /8 Mc /8 / a 105/i fc - 230y7i^7; + 45Mc \ / , 

^ ( 0) ~ ^/ 4 (v^ + v^r { + 47 16^ «^ + v^) J ^ ( }) - 



Radiative decays of the excited heavy baryons might also be measured. Due to the small available phase space, 
radiative decays are expected to have a large branching ratio for excited A c and A b baryons; moreover, they may be 
the dominant decay channels of the excited Af, baryons. It was shown in a model-independent way in Ref. p7| , that 
in the combined limit the electromagnetic decays of the doublet (Aqi, Aq-J in the charm and bottom sectors are 
dominated by dipole radiation. The total decay rate averaged over the initial state and summed over all final states 
is the same for Aqi and Aq 1 [0. The total decay rates in the charm and bottom sector up to NLO are found to be: 

r(A cl ->A c7 ) = l e 2 ^ W f ~™ N ) (I --^=)(1 hC%\),. ,19, 




r(A bl ^A bl ) = ^J ms + mN ) ( 1--^=),! +c>(A)i. ,20, 

where e is the electromagnetic coupling constant (e 2 w 1/137). It is interesting to note that the ratio of these decay 
rates at LO is r(A c i — > A c 7)/r(Afei — > A b j) w 0.2, which significantly differs from 1 — the HQET prediction at 
leading order in 1/ibq. The deviation is due to the large recoil of the heavy quark against the brown muck in the 
charm system where jrtiQ w 1/2. 

At leading order, the values of the heavy baryon observables, eqs (Q), (||), (|l3|), (|l4|), (|T^), (19) and (20), depend 



on one parameter — k. This parameter can be eliminated at LO in a model-independent way by using the mass 
splitting between the ground state and the first orbitally excited state of A c in eq. (0). At LO, the model independent 
predictions are given in Table. [|. The theoretical uncertainty of these model-independent predictions can be large; 
corrections are nominally of relative order A 1 / 2 ~ (1/3) 1 / 2 ~ 60%. To improve accuracy, the next-to-leading order 
corrections must be included. The predictions at NLO contain theoretical uncertainty of relative order A which is 
nominally ~ 1/3 ~ 30 %. 

At NLO, heavy baryon observables, eqs. @j, (||), (|l~5|), (|l4|), (|l^), ( |l9| ) and (20), contain an additional constant 



a. 



Model-independent predictions at this order can, in principle, be made by eliminating k and a in terms of any pair of 
observables. In order to experimentally test these predictions apart from the presently known value of the excitation 
energy of A c (Am c w 330 MeV), another observable (eg. a derivative of the semileptonic decay rate, p or a, or the 
total radiative decay rate of the first excited state of A c ) has to be measured. The sensitivity of different observables 
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TABLE I. Model-independent predictions at LO. The observables are denned in eqs. (Q), (g), (jl|), @, (|l|), ([TJ) and @. 
Corrections are of relative order A 1//2 . 



Fitted 
Observable 


Am b (MeV) 


e 


I 

P (MeV~ 3/2 ) 


D redictions 

a(MeV~ 3/4 ) 


T c (MeV) 


r b (MeV) 


Am c w 330 MeV 


300 


PS 1 


-1.2 x 10 -4 


0.011 


0.025 


0.13 



to the value of a varies. One therefore has to judiciously choose an additional observable. The next-to-leading order 
contribution to such observables should be large enough numerically so as to be stable against the NNLO corrections. 
Based on the sensitivity to a, we will argue that the best observables (in addition to Am c ) are either the second 
derivative of the form factor 0(z) at zero recoil — p — or the slope of the form factor 2(z) at zero recoil — a — or the 
total radiative decay rate of the first excited state of A c or A;, — r c or IV 

To estimate the size of the NLO corrections we can rescale the constants k and a and the observables in eqs. (0), 
(||) , (|lj) , (|l4|) , (|l^) , ( |l9| ) and ( |20| ) using the typical momentum scale of the collective degrees of freedom. It is natural 
to define this scale by, 

A= (fi c Am c 2 ) 1/3 & 410 MeV, (21) 

with A 3 being equal to the value of k determined from the leading order term in eq. (^): n — fi c Am c 2 (l + ©(A 1 / 2 ). 
In terms of the scale A, dimensionless constants R and a can be defined by, 

k = A- 3 k, a = A~ 5 a. (22) 

The dimensionless constant R is unity when only leading order terms in the combined expansion of the observables in 
eqs. (0), (g), dT|), Q, @, @ and @ are kept. When NLO terms are included, the constant R has the form: 

R = (1 + 5R) (1 + 0(A)) , (23) 

where SR is of relative order A 1 / 2 . Hence, when working up to NLO, the LO terms containing R may be linearized with 



respect to SR. Moreover, since in all of the observables in eqs. (0), (||), (|l3|), (|l4|), (|lg|), (19) and fl2C| ) the constants 
k and a contribute only via the ratio a/ n n (where n is an integer or half-integer), the value of R in these expressions 
can be taken to be unity when NNLO corrections are neglected. In addition to using dimensionless constants R and 
a, we can rescale the observables according to their dimensions using appropriate powers of A. The corresponding 
dimensionless observables are given by, 

Am.EAm.A- 1 = (1 + -5R) J — + - } — + 0(A 3 / 2 ) 

2 y fi c 4! p c 

« (0.80 (1 + ~5R) + 0.13 a) (1 + 0(A)) , 

Am b = Am b A" 1 = (1 + ~6r)J — + ^— + 0(A 3 / 2 ) 

i y fib 4\ /it 

« (0.72 (1 + ]-8R) + 0.11 a) (1 + 0(A)) , 

MB,. ( **fj*L fl + ^ ^-/f? ) (1 + OW) 
(v^+VMc) V 4! 16VMfcA 1 c(V^fc + VJtc) J 

w (0.998 + 5.6 x 10" 4 a) (1 + 0(A)) , 

w (-0.998 (1 - -<5k) + 0.15 a) (1 + 0(A)) , 

ff ^3/ 4 _ 4 A*b /8 ^c /8 A 1^_ a 105/x b - 230^6/^c + 45^ c 



rr = rr A 3 ' 4 - ^ ^ 1 - -J)k -4- ^ ^y^c . !r ^c , , 

"A3/4( V ^+^)2^ 1 4 4! lQA 1 / 2 y /JIbJI^(^/]Ib + y/J^c) ) hC ' An 

» (0.969 (1 - -5k) - 0.07 a) (1 + 0(A)) , 
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TABLE II. Model-independent predictions at NLO. Observables are given in eqs. (@), (§), @, @, © and @. 

The typical scale A is denned as A = (/icAm,; 2 ) 1 ' 3 ~ 410 MeV. Corrections are of relative order A. 



Predictions 


Am c and p 


Fitted Observables 

Am c and a 


Am c and V c 


Am b A -1 


1.29 + 0.57 (pA 3/2 ) 


0.77 - 0.05 (o-A 3/4 ) 


0.71 + 0.27 (r c A _1 ) 


e 


0.95 - 0.05 (pA 3/2 ) 


0.99 + 0.004 (aA 3/4 ) 


0.999 - 0.02 (r c A~ 4 ) 


pA 3 ' 2 




-0.92 - 0.08 (<rA 3/4 ) 


-1.01 + 0.46 (r c A^ 4 ) 




0.25 - 0.72 (pA 3/2 ) 




0.98 - 0.33 (r c A^ 4 ) 


TcA- 1 


(2.67 + 2.62 (pA 3/2 )) x 10" 3 


(2.67- 2.15 (aA 3/4 )) x 10~ 4 




r^A- 1 


(1.27 + 1.24 (pA 3/2 )) x 10 -2 


(1.30 - 1.02 (crA 3/4 )) x 10" 3 


(0.029 + 5.77 (rcA^ 1 )) x 10 -3 



f c . %Y C A-i = ( m^L)\* ( 1 + SR _ * + (X)) 

e z \ m B m N ) \ 4! ^/JT C J 

« (0.05 (1 + 8k) - 0.01 a) (1 + 0{\)) , 

f 6 - ^r b A- = f !»£±™ V A a f i + te _ ^ * \ + 

V m s m N J \ 4! ^/Ttfc / 

w (0.26 (1 + - 0.04 a) (1 + 0(A)) , (24) 

where we linearized the terms containing R. In the NLO terms this amounts to taking R to be 1 as was discussed 
above. In the second set of equalities in eq. ( p4| ) the experimental values for the masses have been used. 

The rescaled observables in eq. (|24|) depend linearly on the dimensionless constants 5R and a. The value of a is 
expected to be of order unity. The NLO corrections for all observables in eq. (pi|), except for Oo, are 7% to 20% 
when a is taken to be unity. The NLO contributions to the the value of Oo is suppressed by four orders of magnitude. 
The leading order value of ©o is very close to unity (see Table |) — the leading order value predicted by HQET. 
The deviation from unity is an exceptionally poor way to extract a since the sensitivity is so small numerically that 
one expects higher order terms to be significant. The sensitivities with respect to R can be used to decide which 
observables are best suited to extract values of R and a. 

Formally the 1/toq and 1/N C corrections in the combined expansion are of the same order. However, the heavy 
quark expansion is generically expected to be more a accurate approximation for heavy baryons than the large N c 
expansion. If Amf, is used apart from Am c to eliminate k and a, the resulting relations contain mass differences, 
Am c — Ami,. In the heavy quark limit this difference is zero up to 1/toq corrections. At leading order in the combined 
expansion this difference is approximately 30 MeV (Table ||) . The NNLO corrections which are beyond the accuracy 
we are working to can be estimated from the mass splitting of the states of the doublet (A cl , A*i). This splitting — 
approximately 40 MeV — is the NNLO effect containing 1/mq and 1/N C corrections. Hence, predictions containing 
Am c — Arrib and Am c + Anib are sensitive to NNLO corrections and should not be used to extract k and a. 

Any of the remaining observables — p, a and r c (or Tb) — can be used along with Am c to eliminate k and a. 
Using Am c and p as primary observables, we can obtain model-independent predictions for other observables. These 
predictions are listed in the second column in Table [Fl|. Similarly, using eq. (p4]), R and a can be extracted from pairs 
(Am c , a) and (Am c , T c ). Corresponding model-independent relations are given in the third and the fourth columns 
in Table [n]. The corrections to relations in Table [n] are generically of order A which nominally is 1/3 ~ 30%. In 
order to test NLO predictions, measurements of one of the observables — the curvature of the form factor Q(z) at zero 
recoil, the slope of the form factor 3(z) at zero recoil, or the total radiative decay rate of the first excited state of 
A c — has to be measured. 

We have discussed the model-independent relations between observables which arise in the combined heavy quark 
and large N c expansion at leading and next-to-leading order. These observables were derived in the framework of 
an effective theory in the combined heavy quark and large N c limit in |l(|[l7]]. At leading order in the combined 
expansion a model independent expressions for a number of heavy baryon observables contain one constant. This 
coefficient can be eliminated using the known value of the excitation energy of the first excited state of A c , Am c . The 
resulting predictions up to corrections of order A 1 / 2 are listed in Table |. To obtain predictions at next-to-leading 
order in the combined expansion, an additional observable is required. At the present time the data is not available. 
Observables in eqs. (Q), (g), (p^|), (|l4|), (|l§|), ( |l9| ) have different sensitivities with respect to the next-to- next-to-leading 
order corrections. We have discussed observables (when data becomes available) should be used in addition to the 
excitation energy Am c to predict the remaining observables. The model-independent relations at NLO are given in 
Table [n[ 
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